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1. Introduction 

Quantum Chromodynamics (QCD) is the "nearly perfect" theory of the strong interactions [1]. 
The QCD Lagrangean, expressed in terms of fundamental quark and gluon fields, is rich in sym- 
metries and is parameter free if quarks are taken to be massless. Phenomena in the world where 
the strong interaction holds sway arise as a consequence of the dynamics of the quark and gluon 
fields. In particular, the QCD vacuum plays an important role in the structure of strongly interact- 
ing matter. A subtle example of the role of the vacuum is chiral symmetry breaking, whereby the 
masses of light quarks are leveraged via the Nambu-Goldstone mechanism into the much heavier 
masses of the spectrum of observed hadrons. 

Because QCD is a theory where most phenomena are "emergent" phenomena, it is especially 
interesting to ask what happens to QCD matter when squeezed to ultra-high densities or heated to 
ultra-high temperatures. The role of the vacuum is greatly modified and strongly interacting matter 
may behave very differently under extreme conditions relative to those found at low temperatures 
and densities. 

The structure of QCD has immediate consequences for extremely hot and dense matter which 
resolve puzzles about the strong interactions present since its early days such as, for example, "lim- 
iting temperatures". Because of asymptotic freedom, QCD matter asymptotically must be a decon- 
fined system of weakly interacting quarks and gluons [2]. Likewise, for matter at extremely large 
baryon densities, the close packing of hadrons gives way to quark matter [3]. Because screening of 
quarks and gluons plays a significant role in the dynamics of QCD matter at high temperatures and 
densities, such a deconfined state is called a Quark-Gluon Plasma [4]. Our universe was a QGP 
~ 10~ 5 seconds after the Big Bang. 

Another key feature of QCD is "infrared slavery" at large separations (or low temperatures and 
densities)-a phenomenon closely related to confinement. As a consequence, the QGP undergoes 
a transition from deconfined quark-gluon matter to confined hadronic matter. Similarly, chiral 
symmetry which is restored at high temperatures is broken at lower temperatures. It is therefore 
clear that the theory must contain a rich phase diagram for variations in temperature and baryon 
chemical potential. Significant theoretical developments in lattice gauge theory and an explosion of 
computing power have made it possible to explore aspects of the phase diagram in great detail. We 
refer the reader to Owe Philipsen's talk at this conference for details on the latest developments [5]. 

Our focus here is on heavy ion collisions at ultrarelativistic energies where the aim is to pro- 
duce deconfined QCD matter in the laboratory-this matter is the hottest and densest matter pro- 
duced on earth and lasts for ~ 3 • 10~ 23 seconds. We will argue here that a wide range of results 
from the RHIC experiments coupled with theoretical developments have given rise to a picture of 
heavy ion collisions on which there is a sufficiently broad consensus for us to describe it as a "stan- 
dard model" of heavy ion collisions. Formulating a standard model is useful and perhaps essential 
for an intellectually coherent narrative of the science, to expose clearly possible weak chinks in 
such armor and to bring clarity to future paths along which the field could develop. 

2. A standard model of heavy ion collisions 

Exploring the QCD phase diagram at high temperatures and densities in heavy ion collisions 
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is challenging because the collisions are intensely dynamical processes where contributions to any 
measured final state can arise from different regions in the space-time evolution of the collision. 
One therefore needs to understand the entire chronology of the collision starting from the properties 
of the high energy nuclear wavef unctions. While historically the focus of heavy ion collisions was 
primarily on understanding the thermal properties of deconfined QCD matter, there has been a 
gradual realization that other aspects of the collision open additional interesting windows into the 
many-body dynamics of QCD in collision. 

II ! • 

CGC initial Glasma 

Singularity 

Figure 1: A standard model of the space-time evolution of heavy ion collisions. 

A standard heavy ion model consists of the chronology shown in fig. 1 . Before the collision, 
the nuclear wavefunctions are described by correlated multi-parton states called Color Glass Con- 
densates (CGC). The CGCs shatter in the collision producing deconfined non-equilibrium QCD 
matter called the Glasma. This Glasma likely thermalizes to form a strongly correlated Quark- 
Gluon Plasma (sQGP), which subsequently undergoes a cross-over transition to hadronic matter 
with much lower energy densities. The sQGP flows like a nearly perfect fluid as measured by a 
very low value for the dimensionless ratio of the shear viscosity to entropy density in the fluid. 
When the expansion rate of the system exceeds the scattering rate of the hadronic constituents, the 
matter free streams to the detectors. It is sometimes opined that the earliest stages of the nuclear 
collision are less likely to be understood than the latter stages. My view is that precisely the op- 
posite may be true-one has some confidence that weak coupling techniques in QCD can be used 
to systematically compute features of early time dynamics. As the system evolves, weak coupling 
is less reliable and systematic computations can perhaps be performed for a limited sub-set of the 
system's properties 1 . Nevertheless, as we shall discuss later, a number of theoretical techniques, in 
combination with experimental results, can shed considerable light on many-body dynamics in the 
sQGP. 

In the following, we shall briefly summarize the theoretical ideas and the extant empirical 
evidence for each of the stages of the standard heavy ion model. It is impossible to do full justice 
to this topic in the limited space available for this talk and I will therefore have to frequently refer 
the reader to the available literature. 

3. Before the collision: wee parton correlations in the Color Glass Condensate 

The appropriate asymptotics for multi-particle production in QCD is the Regge-Gribov asymp- 

'This does not mean we have all the answers for early time dynamics in weak coupling. An example is thermaliza- 
tion which is still poorly understood. However, it is less well understood in strong coupling frameworks. 
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totics of the theory corresponding to Q 2 = fixed, x^ — > and s — > °°. In the framework of the theory 
where the parton model arises, this limit corresponds to very high partem phase space densities. 
In contrast, the more highly developed Bjorken asymptotics of the theory, Q 2 — > °°, s — > °° with 
XBj = fixed, corresponds to low occupation numbers and is more appropriate for the description of 
rare high momentum transfer processes. A natural consequence of Regge-Gribov asymptotics at 
small x is parton saturation [6], which, in the light front formulation of QCD, corresponds to maxi- 
mal occupancy of ~ I /as for gauge field modes below a momentum scale Qs, the saturation scale. 
An equivalent description [7] of parton saturation in the hadron rest frame corresponds to a small 
size (dipole) probe of size ~ 2/Qs scattering off a hotspot within the hadron with unit probability. 

The high occupancy of saturated gluon modes (for Q 2 S 3> Aqcd) suggests that the nuclear 
wavefunction at high energies (or small x) can be simply described by classical fields [8]. Further- 
more, because Qs is large, this high occupancy state is weakly coupled and its very non-perturbative 
behavior can be computed systematically order-by-order in as beyond the leading 1 /as classical 
contribution. The Wilsonian renormalization group (RG) equations that describe the evolution of 
this high occupancy state with energy are stochastic in nature, which explains why it is described 
to be a Color Glass Condensate (CGC) [9]. These RG equations, called JIMWLK equations [10], 
govern the energy evolution of n-point wee parton correlations in the hadron wavefunction. These 
carry much more information about QCD dynamics than parton distribution functions and are in- 
teresting to study in their own right. 

What is the evidence for parton saturation and the CGC ? An early remarkable phenomenon 
was geometrical scaling of the inclusive cross-section with Q 2 /Q s [1 1]- Phenomenological models 
that incorporate the physics of parton saturation do a very good job of describing HERA inclusive, 
diffractive and exclusive final states [12], fixed target e+A inclusive data [13] and RHIC A+A 
multiplicities [14] with an economical and consistent set of parameters. An early success of the 
saturation picture at RHIC was the description of the centrality dependence of the multiplicity 
distribution in A+A collisions [15]. The phenomenological saturation models have been shown to 
give a good description of the LHC p+p data including the n-particle multiplicity distribution[16, 
17]. Fig. 2 (left) shows the saturation scale versus impact parameter for three different values of x 
in p+p collisions extracted from HERA data in two different saturation models. The figure on the 
right shows the saturation scale as a function of x and A. Both plots show the saturation scale in 
the fundamental representation. In the adjoint, relevant for hadronic collisions, one multiplies the 
value of Q 2 S shown by 9/4. 

While the phenomenological models work very well, it is more satisfactory to confront next-to- 
leading order CGC predictions with the data. There has been quite a bit of progress in this direction. 
The Balitsky-Kovchegov (BK) [18] equation, the large N c "mean field" equation in the CGC for 
the dipole forward scattering amplitude (corresponding to a 2-point Wilson line correlator) has now 
been computed at the next- to-leading log (NLL) level [19]. An evolution equation including the 
running coupling corrections in the NLL expressions [20] gives good agreement with the HERA 
inclusive data [21], e+A inclusive data [49] and the forward single inclusive deuteron+gold data 
from RHIC [22]. More interestingly, it describes [23] the systematics of the disappearance of the 
away side hadron in di-pion production in deuteron-gold collisions at forward rapidities observed 
by the STAR collaboration [24]. While theoretical caveats [25, 26] to the predicted result may be 
quantitatively important, they are unlikely to alter its qualitative features. A more serious obstacle 
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Figure 2: Left: The saturation scale in the fundamental representation as a function of impact parameter 
extracted from HERA data in two different saturation models. Right: Picture of hadrons at high energies as 
resolved by a lepton probe. 



to this interpretation is from the possibility [27] that a azimuthal angle independent pedestal effect 
could mask this result. A clean experiment which will be able to settle the issue conclusively is to 
measure the same final state at a future electron-ion collider [28]. 

In A+A collisions, most saturation models made predictions for the LHC that were on the low 
side by ~ 20% [29] compared to the recently released ALICE data [30] for Pb+Pb collisions at 2.76 
TeV/nucleon. Integrated multiplicities themselves are not terribly sensitive to dynamics. A 20% 
discrepancy in the multiplicity corresponds to a corresponding discrepancy of 10% in Q$, which 
can easily be accounted for in revised fits. Indeed, when nuclear geometry effects are included, 
the NLO-BK model predicted a centrality dependence of the multiplicity [31] that shows excellent 
agreement with the ALICE data. 

While each of these comparisons of CGC/saturation models to data in e+p/A, p+p, d+A and 
A+A collisions can be (and are) challenged, taken together they provide a remarkably consistent 
and detailed picture of a significant body of data suggesting that saturation effects have been seen 
and are essential to understanding multi-particle production at high energy colliders. At the LHC 
(and in forward observables at RHIC), the saturation scale is large enough that weak coupling 
techniques are applicable. This therefore suggests the exciting possibility that physics previously 
thought inaccessible is now amenable to systematic computation. 

4. The little bang and the Glasma 

It was understood early on in QCD that multi-particle production in hadron-hadron collisions 
arises from the scattering of their wee parton clouds [32]. In the CGC, these are classical fields with 
their field strengths localized on the Lorentz contracted widths of the nuclei, so to lowest order in 
the coupling constant (of order I /as), the scattering of two CGCs is described by the collision of 
two classical fields. This is equivalent to solving Yang-Mills equations with initial conditions deter- 
mined by the classical fields for each nucleus [33]. These classical equations have been solved and 
the single inclusive [34] and more recently the double inclusive [35] distributions determined. To 
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this lowest order, the classical equations (and particle distributions) are boost invariant. The gauge 
field configurations are color screened on distances of order l/Qs and right after the collisions 
are dominated by longitudinal color-electric and color-magnetic fields that carrry Chern-Simons 
charge [36]. The QCD matter corresponding to early time non-equilibrium dynamics has been 
dubbed "the Glasma" [37]. 

The classical configurations in the Glasma are far from the whole story. Quantum fluctuations 
profoundly enrich this picture 2 . Before the collision, quantum fluctuations generate large loga- 
rithms as\&{\/x\fl) that are °f the same order as the leading classical term and therefore have to 
be resummed at each order in perturbation theory. At each loop order, there are also large multiple 
scattering contributions that need to be summed to all orders in dsn, where n ~ l/o^ is the den- 
sity of parton sources. High energy factorization theorems proven recently for inclusive quantities 
show that the leading logs in x and multiple scattering contributions can be factorized 3 from the 
final state dynamics into universal density matrices that satisfy the JIMWLK equation [39] con- 
voluted with the inclusive observable computed at leading order. Thus correlators of Wilson lines 
extracted from e+A scattering can in principle be used as inputs into computing A+A final states. 

The high energy factorization theorems are a powerful tool to compute the early stage space- 
time evolution of quantities such as the stress-energy tensor or correlators of the stress-energy 
tensor ab initio in heavy ion collisions. They allow one to compute for instance long range rapid- 
ity correlations of inclusive observables that are sensitive to the earliest times in the collision. A 
consequence of the formalism is the Glasma flux tube picture [40] where n-particle correlations 
are simply proportional to {S ^ / {I / Q 2 s )) l ^ n for n>2, the ratio of the nuclear transverse area S± 
to the area of a flux tube 1 / Q 2 S . These geometrical correlations give rise [41] to the negative bino- 
mial distribution which explains multiplicity distributions [17] and forward-backward multiplicity 
correlations [42] in hadronic collisions. 

Long range rapidity correlations that are sharply collimated around AO « were observed at 
RHIC [43]. These structures, called "ridges" based on their visual appearance, are present both 
for triggered and untriggered two particle correlations. These ridges are not present in peripheral 
events and are most prominent in the most central events. Long range rapidity correlations are in- 
teresting because they are a "chronometer" of strong color field dynamics at early times in hadronic 
collisions. This is illustrated in fig. 3 (left). 

In the Glasma flux tube picture [40], the ridge phenomenon is a consequence of the long range 
rapidity correlations in flux tubes of transverse size l/Qs formed early in the collision that are 
subsequently boosted in the final state by radial flow [45]. The radial flow provides the near side 
angular collimation. The combination of flux tube structures in the initial state and radial flow 
gives an excellent description of the published RHIC data [46, 47]. The ALICE collaboration has 
presented preliminary evidence of a ridge in Pb+Pb collisions [48]-we would expect the amplitude 
to follow the trend seen at RHIC and be larger than the RHIC values. At the LHC, one expects the 
ridge correlations to show deviations from boost invariance [49] for AT7 > 4 units and the amplitude 
of the ridge to decrease with increasing trans vere momentum for the pairs [53]. 

2 This is fortuitous because the naive leading order picture is not what is seen in experiments. 
3 The underlying basis for this factorization is the observation that the computation of inclusive quantities in field 
theories with strong time dependent sources can be formulated as an initial value problem [38], 
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Figure 3: Left: Plot illustrates how long range rapidity correlations are causally limited to provide infor- 
mation on early time dynamics. Right: A collimated structure around A<p s» that is long range in rapidity 
seen by the CMS collaboration in high multiplicity p+p collisions at c.m.s energies of 7 TeV. This ridge-like 
structure is not seen in Monte-Carlo models while the away side peak at A(j) ss n is well reproduced in these 
models. 

In a recent paper, the CMS collaboration announced the discovery of a ridge in high multiplic- 
ity events in p+p collisions [54] shown in fig. 3 (right). This effect was predicted [50] based on our 
formalism [49] but was not advertised because it was a small effect. We subsequently showed the 
Glasma flux tubes predictions are in qualitative [51] and quantitative agreement with the experi- 
mental data [52]. We note that long range forward-backward correlations were predicted for p+p 
collisions at the LHC in a related framework [55]. 

Another effect that is very sensitive to properties of the Glasma is the Chiral Magnetic Effect 
(CME) [56]. This is clever idea that suggests that sphaleron transitions in deconfined matter subject 
to an external magnetic field can lead to charge separation in the direction of the magnetic field that 
locally breaks P and CP though of course it is preserved in the event as a whole. The STAR collab- 
oration published data on same and opposite charge separations which appeared consistent with the 
CME [57]. Several recent papers however suggest alternative analyses of the STAR data [58]-the 
RHIC low energy run may help clarify which interpretation is correct. Nevertheless, the CME has 
stirred a lot of theoretical interest and has ramifications for fields outside heavy ion collisions [59]. 

How the Glasma thermalizes to a QGP is an outstanding problem that has not been solved even 
in weak coupling despite much theoretical progress in recent years. An important ingredient is the 
role of instabilities [60]. In the CGC framework, unstable quantum fluctuations play an essential 
role in close analogy to the situation in inflationary cosmology [61]. Rapidly growing quantum 
fluctuations can be resummed; in a toy scalar theory, the corresponding energy-momentum tensor 
obeys ideal relativistic hydrodynamics [62]. If these considerations can be extended to a gauge 
theory, one could have "hydrodynamic-like" flow without early thermalization. 

5. The perfect fluid 

A key result of the RHIC experiments is the large flow measured suggesting that hydrodynam- 
ics is applicable. Hydrodynamics is the right effective field theory to describe the long wavelength, 
late time behavior of quantum field theory. The surprise at RHIC is that it works at much earlier 
times than simple estimates would suggest. A powerful measure of the degree of flow is v%, the sec- 
ond moment of the anisotropy in the single particle distribution. Specifically, it measures how effi- 
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ciently hot matter converts spatial anisotropics in the initial distribution to momentum anisotropics 
at later times-the most efficient way to do this is by applying ideal relativistic hydrodynamics. 

A quantitative measure of the efficiency of flow is given by the ratio T7 /s of the shear viscosity 
to the entropy density in the fluid. For simplified boost-invariant ("Bjorken") hydrodynamics, the 
viscous contribution to the evolution in proper time % of the energy density relative to the ideal 
term is given by r\/s/xT\ because XjxT ~ 0(1), this contribution is small when r\/s << 1. In 
kinetic theory, r\/s « T r ei a x.Aq Uant ., i n units of h/ks, where T re iax. is the momentum relaxation time 
in a fluid in response to a shear stress and T qua nt. is the thermal Compton wavelength in the fluid. 
One therefore expects naively a lower bound of unity on 77 js. However, this expectation is not 
correct because there is no reason kinetic theory is relevant in this domain. In fact, for strongy 
correlated systems one expects kinetic theory to fail 4 . Nevertheless, on very general grounds one 
expects a lower bound on r\ /s. Is there such a bound and how does the value extracted from RHIC 
experiments compare ? 

While the answer to the first part of the question is not known definitively, there is a conjectured 
lower bound 5 of r]/s = I /An [63]. This bound was derived for N=4 supersymmetric Yang-Mills 
(SYM) theory and makes use of the Maldacena AdS/CFT conjecture relating SYM field theory in 4 
dimensions to 10 dimensional classical gravity in a background of D3 branes [65]. Specifically, the 
viscosity in SYM is related to the absorption cross-section of a graviton on a black brane, which 
like the entropy density, is proportional to the area of the brane. 

What is T}/s extracted empirically ? For this one relies on viscous hydrodynamic models-the 
algorithms for these have developed significantly [66]. While there are uncertainties related to 
the initial conditions and other details of the numerical implementation, its probably safe to say 
that rj/s < 0.4 [67]. A comparison of r]/s in various fluids at temperatures close to the critical 
temperature is shown in fig. 4 (left). These values of r]/s are much lower than that of water or even 
liquid helium but are comparable to those of strongly correlated lithium-6 atoms in the unitarity 
limit. It is remarkable that two systems whose temperatures differ by 10 19 and viscosities by 10 26 
flow nearly identically. 

Most weak coupling estimates of tj/s in the QGP give values for 77 js that are significantly 
larger [68]. Does this rule out a weak coupling description of flow completely? One way out is if 
the system had a small "anomalous" viscosity (as seen in electromagnetic plasmas) which mimics 
a small kinematic shear viscosity [71]. If, as discussed previously, quantum corrections to the 
Glasma lead to ideal hydrodynamic behavior, this could provide a quantitative mechanism for the 
"anomalous" early time viscosity. Detailed studies are feasible in near future on how much the 
empirical lower bound on 77/5 can be moved upwards. 

The ALICE collaboration has released first results on vi seen in Pb+Pb collisions at ^]~s = 
2.76 TeV/nucleon [72]. The integrated V2 is about 30% higher than at RHIC but V2(p±) is nearly 
identical to that seen at RHIC. This is what one would expect if a nearly perfect fluid had been 
formed at RHIC and is consistent with the predictions from viscous hydro models [73]-see fig. 4 
(right). If significant entropy were generated, it should impact the centrality dependence of the 

4 Albeit, one should note that kinetic theory often does much better (as is the case for cold atomic gases) than its 
condition of applicability would suggest. 

5 There are recent suggestions that in certain higher derivative gravity theories an even lower bound may be at- 
tained [64]. 
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Figure 4: Left: Behavior of rj /s in a number of fluids measured relative to deviation from the critical 
temperature [69]. Right: Data on v% from the ALICE collaboration [70] compared to predictions from 
viscous hydrodynamic models. 

multiplicity distribution, which is consistent with initial state predictions for both RHIC and the 
LHC. This centrality dependence also has the potential of constraining the thermalization time of 
the Glasma [74]. Interesting measurements at RHIC (and in near future at LHC) on V2 fluctuations 
and higher moments of the anisotropy distribution are quite sensitive to fine details of the initial 
conditions and final state effects [75]. For instance, the ratio v^/v^ gives very different results for 
different ansatze of the viscous corrections to the single particle distributions at freeze-out [76]. 

6. Hard probes of QCD matter 

Hydrodynamics paints broad brush strokes of the dynamics of strongly correlated matter. To 
obtain further insight, one looks to hard probes of the QCD medium. Understanding the dynamics 
of the medium from these is very interesting albeit very challenging as well because a colored hard 
probe is typically sensitive to the entire time history of the system. At RHIC, it was observed that 
7?,u-the ratio (normalized by the number of binary nucleon-nucleon collisions) of the inclusive 
hadron spectrum in A+A collisions relative to the same quantity in p+p collisions-was suppressed 
by a factor of 5 for % s and is nearly flat out to the highest /jj_'s measured at RHIC [77]. A control 
deuteron+gold measurement [78] established that this phenomenon was a final state effect arising 
from the interaction of partons with the colored medium. 

"Jet quenching" was suggested as a probe of the QGP a long time ago [79] and radiative 
energy loss by partons traversing the medium was identified as a weak coupling mechanism that 
would explain it [80, 81]. The energy loss of partons can be related to a transport coefficient q 
that characterizes the properties of the medium. There are several formalisms based on the original 
works that have been developed to confront the heavy ion data on jet quenching [82]. These differ 
quite a bit in detail and give values of q that differ considerably-for a clear and interesting recent 
discussion of the sources of some of these differences, see Ref. [83]. 

The flat behavior of Raa, the large values of V2(p±) out to large p±, and the unexpectedly 
large energy loss of heavy quarks are difficult to reconcile in a weak coupling framework without 
fine tuning. An example of a fine tuning argument is one that explains the flat behavior of Raa at 
RHIC as arising from a combination of effects, in particular phase space constraints at large xj, 



9 



Standard model of heavy ion collisions 



Raju Venugopalan 




Figure 5: Left: Data from the ATLAS collaboration showing significant jet asymmetry for central Pb+Pb 
collisions. The data also shows that despite the asymmetry, the jets are back-to-back. Right: Next-to-leading 
order estimate compared to leading order estimate of contributions to the second moment (q) of the elastic 
scattering rate of a hard jet with a hot medium for transverse momenta close to the medium temperature. 

that modify the weak coupling prediction that the suppression goes away with increasing p±. This 
should be different at the LHC where there is no large xj constraint. The first Raa data at LHC 
for charged hadrons indeed appear to favor this interpretation [84]. One would like however to 
see the spectrum of 7T°'s and a wider range in p± before more definitive statements are made. An 
important factor that should be kept in mind while interpreting Raa is the very large initial state 
gluon shadowing at LHC energies [85, 22]. 

The apparent strong opacity of the QGP has triggered a considerable amount of work [86] 
to explain light and heavy quark energy loss in the same AdS/CFT framework that gave us the 
lower bound on v\ /s. There are two sorts of treatments here. In one, the jet production mecha- 
nism is treated perturbatively and factorized from the non-perturbative interaction of the jet with 
the medium [87]. The latter is expressed in terms of correlators of Wilson lines which is com- 
puted using the AdS/CFT correspondence. In the other, jet quenching is treated completely in the 
AdS/CFT framework, typically by looking at the energy loss of heavy quarks [88]. QCD is how- 
ever significantly different from N=4 SYM, especially in the region around T c , where one notices 
strong conformal symmetry breaking effects. Until one has a good candidate for a QCD dual in the 
AdS/CFT framework, perhaps its value is primarily to provide guidance into dynamical questions 
that one cannot easily answer in the QCD plasma in the region where weak coupling is suspect. 
For other approaches in addressing this question, see for instance ref. [89]. 

It is important to understand fully the limits of the applicability of the weak coupling frame- 
work. For instance, present treatments of the partem interaction with the medium are at tree level. 
An interesting quantity to compute is the elastic scattering rate for a hard particle in the medium; 
the transport coefficient q is the second moment of this quantity. The first NLO 0(g) computation 
of this quantity has been reported recently [90] ; the correction to the leading order result is signif- 
icant. The estimate is shown in fig. 5 (right). This raises the question whether resummations a la 
the pressure in finite temperature field theory are required. 

As promised, jets have been produced in copious amounts in Pb+Pb collisions at the LHC. The 
ATLAS and CMS collaborations confirm observation of jet quenching in di-jets [91, 92]. Looking 
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at trigger jets with Ej\ > 100 GeV and away side jets in the opposite hemisphere with £"7-2 > 25 
GeV as a function of the variable Aj = {Ej\ — Er%) /Ej\+ Et2, they find a centrality dependent di- 
jet asymmetry that is strongest for the most central collisions. At the same time, the vast majority 
of the di-jets are back-to-back-see fig. 5 (left). A simple explanation [93] that makes intuitive sense 
follows from the observation that firstly, the softest components of the jet shower decohere sooner 
and secondly, are more easily deflected out of the jet cone by scattering off the hot medium thereby 
degrading the energy of the jet, while preserving its angular structure. The first estimates of q are 
consistent with a weak coupling scenario. Clearly, a lot more data can be anticipated in the very 
near future that will add considerable depth to the existing picture of jet quenching. 

7. Looking ahead 

I have outlined here the bare bones of a standard model of heavy ion collisions. Even within 
this structure, there is much that we don't understand, but we have a better idea of what that is. 
With LHC very quickly showing us what is possible and the high luminosity and detector upgrades 
at RHIC on track, we can be optimistic the outstanding issues will be solved. Because many-body 
QCD is a rich and subtle science, a natural extension of on-going research on hot matter is to 
explore the complementary many body dynamics of "cold" nuclear matter with a high energy, high 
luminosity electron-ion collider. 
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